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Abstract 

In this paper, we introduce the SD-Quantum Stationary Hamilton Jacobi 
Equation for a central potential, and establish the 3D quantum law of 
motion of an electron in the presence of such a potential. We establish 
a system of three differential equations from which, and as a numerical 
application, we plot the 3D quantum trajectories of Hydrogen's electron. 
We show also the existence of nodes for these trajectories. 
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1. Introduction 



Three years ago, we presented with Bouda a deterministic approach of quan- 
tum mechanics by which we attempted to restore an objective reahty of quantum 
phenomena. This approach consist on the introduction of a ID quantum La- 
grangian from which we derived the ID quantum law of motion Then, we 
plotted the quantum trajectories for several one dimensional systems After, 
we generalized our approach to the relativistic systems [3], and obtained the 
relativistic quantum trajectories and established the existence of nodes also for 
such systems 0. we also presented a study of the relativistic spinning particle 
using the deterministic approach that we developed in Refs. PEIEIEI- Indeed, 
we gave, first, the relativistic quantum stationary Hamilton-Jacobi equation 
for a particle with a half (1/2) spin RQSHJESi then, we established the 
Quantum law of motion of such a system |^. 

Our deterministic quantum approach is interesting, however it was elaborated 
only for one dimensional systems. This fact was a real handicap for studying 
the nature since the most phenomena happens in real 3D space. For this reason, 
we presented in Ref. [7j , a generalization of our approach to three dimensional 
systems. First, we set the expression 

So(r) = ^larctan (a + b] (1) 

V nr) J 

of the quantum reduced action as a solution of the 3D-QSHJE 

V. (i?2(f)v5o) - . (6) (2) 

The expression of the function R{r) is given as 

[{ae + bcj,f + ^^]-^ . (3) 

Note that 6 and (j) are two real independent solutions of the 3D Schrodinger 
equation 

A* + y(r)* ^ E^. (4) 

zniQ 

a, b and k are real constants. After the establishment of the 3D quantum 
reduced action Eq. (1), we introduced a coordinates transformation with which 
the 3D-QSHJE (Eq. (2)) takes the form of the 3D Classical Hamilton-Jacobi 
equation (3D-CSHJE) j7] exactly as what Faraggi and Matone done for the ID- 
QSHJE 8 . Then, we introduced the 3D quantum Lagragian from which we 
derived the relation between the conjugate momentum and the velocity of the 
particle expressed as [7j 

V . VSo ^2[E- V{r)] . (5) 

We consider this last equation as the Quantum law of 3D motions. It is because 
of -as we will see later- the fact that we should use this equation in order to 
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plot the 3D quantum trajectories. Also, it can be seen as a generalization into 
3D of the ID quantum law of motion already established in Ref. 

x^^2{E~V{x)). (6) 
ax 

The elaboration of the 3D quantum law of motion (Eq. (5)) for a general form 
of the potential is an important step to investigate the 3D motions. However, 
we should plot the 3D quantum trajectories when it is possible in order to 
understand the behaviour of the quantum particle in 3D spaces. With this aim, 
we should act exactly as we have done for ID quantum and relativistic quantum 
trajectories |21E1- Indeed, we must first resolve the Schrodinger equation (Eq. 
(4)) to find the functions and (/>, which we replace into the expression of the 
reduced action (Eq. (1)). Then, we replace the expression (2) into Eq. (5). 
After resolving Eq. (5), we can plot the quantum trajectories. Note that for 
a general potential, plotting the 3D quantum trajectories is not evident. This 
is why, in this article, we study the case of the central potential, and plot the 
trajectories of an electron moving under the atomic potential. 

The paper is organized as follows: In Sec. 2, we introduce the 3D-QSIIJE and 
its solutions for the central potential. Then, in Sec. 3, we establish the quantum 
law of motion for this potential and elaborate a system of three differential 
equations which describe the dynamical behaviour of the quantum particles. In 
Sec. 4, we plot the 3D quantum trajectories of the Hydrogen's electron, for 
several bound states. Finally, In sec. 5, we present a conclusion in which we 
discuss our results. 

Some calculation and demonstrations are given in the and of the article as 
appendix. 

2. The 3D-QSHJE for the central potential 

The central potential is a potential for which only the radial position deter- 
mine its value 

V{r) = V{r) . 

In order to construct the 3D-QSHJE for such a potential, we write the Schro- 
dinger equation (Eq. (4)) with the spherical coordinates (r, i?, Lp) 

*(r, z?, if) is the wave function of the quantum system taken into consideration. 
To resolve Eq. (7), the variables must be separated into 'i/{r,'d,Lp) as [5] 

vl/(r ,^,V5)=i?(r)r(t?)F(^) , 

where R{r), T{'d) and F{ip) are real functions of the variables r, d and (p respec- 
tively. R{r) is the radial wave function, T{§) and F{lp) are the angular wave 
functions. It follows that Eq. (7) will be separated into three equations, each 
corresponds to one variable, 

r'^cPR 2rdR 2mar'^^^ 
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+ m?F = O. (10) 

A is a constant such as 

X = l{l + 1), 

I being a positive integer or vanish, and mi is an integer satisfying to 

—l<mi<l. 

is the eingen value of L^, L representing the angular momentum operator. 
mi represents the eingen value of L^, the operator projection of L along the z 
axis. 

Now, let us write Eqs. (8), (9) and (10) in a form analogous to the one 
dimensional Schrodinger equation. First, we take up Eq. (8) and write 

Rir) = ^ , (11) 

where X{r) is a function of the radius. Then by taking twice the derivative of 
Eq. (11) with respect to r, we obtain 

dR _ldX_ _X_ ^^2^ 
dr r dr 

d^R_ld^X ^dX/dr ^^X 
dr"^ r dr"^ 

substituting Eqs. (12) and (13) into Eq. (8), we get 



-h"^ d'^X L . Xh^ ' 



2too dr^ 



+ 



V{r) + 



2mor^ 



X = EX . (14) 



This last equation has the same form as the one dimensional Schrodinger equa- 
tion for X{r) with an energy E and a fictitious potential 

V (r) = V(r) + ^ . (15) 

Likewise, Eq. (9) can be reduced to the form of the ID Schrodinger equation. 
This can be done by introducing the function T{'&) defined as 

T(a?) = sin5 ^?T(t?) . (16) 
Taking twice the derivative of this last equation with respect to i?, we deduce 
dT _ 1 dT 1 cosi} ^ 

_ 1 rf^T cost9 dT 3 cos^ § 1 T 
d^'' ~ sini a? dd'' sini ^dd^ ^sm^ d ^ 2sini ' 
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If we replace Eqs. (17) and (18) into Eq. (19) we obtain 



d'^T (1/4- m^)^ „ 



(19) 



Remark that Eq. (19) has the form of the one dimensional Schrodinger equation 
with a fictitious potential 



(mf-l/4) 
2too sin^ 1? 



and an energy 



1 h"^ 



(20) 



(21) 



4 2mo 

Remark also that Eq. (10) has the same form as the ID Schrodinger equation 
with a vanishing potential and an energy equal to ?i^/2m. 
Because Eqs. (8), (9) and (10) come to the form of the ID Schrodinger equation. 
So, in order to obtain the 3D-QSHJE, let us write the functions X{r), T(t?) and 
F{ip) as 



X{r) = A{r) 

F{ip) = T]{ip) 



Z(r) 



f3e- 



\:Z{r) 



76 



+ se' 



(22) 
(23) 
(24) 



By replacing Eqs. (22), (23) and (24) respectively in Eqs. (8), (9) and (10), we 
get 



1 



dZ 



2mo V dr 



1 (PA , XTP 

+V{r) + 



2mo A dr^ 



2mo 



E 



d_ 
dr 



dZ_ 
dr 



fdL 



Id^ 



sin^i? 



dd 



dL\ 



3:^U^3T =0, 



dM 
dip 



d-d J 



mfn'^ = 



d_ 

dip 



dM\ _ 
dip J 



= . 



(25) 
(26) 
(27) 
(28) 
(29) 
(30) 



Replacing Eq. (26) into Eq. (25), Eq. (28) into Eq. (27) and Eq. (30) into Eq. 
(29), we find 



1 



dz\- 



- - 



2mo \dr J 4mo 



dL 



{Z,r} + V{r) + 



2mor^ 



= E 



■ 2 , ' n' = {x + -)n' 

sm •& 4 



(31) 
(32) 



5 



[dip] 



■{M,<p} 



(33) 



Eqs. (31), (32) and (33) represent the components of the QSHJE in 3-D for 
a central potential. These equations contain the Schwarzian derivatives of the 
functions Z, L and M with respect to r, •& and ip respectively. Note that, 
by deducing the vahic of A from (31) and the value of mi from (33), then 
substituting in Eq. (32), one obtain 



1 



2too 



dZ\' 



1 



dL 
M 



1 



4mo 



1 



sin^ 'd \ d(p 
1 



{Z,r} + ^ {L,^} + 



fdM 
dip 

{M,ip} 



+ 



V{r) 



8mo 8mo sin^ i? 



E . (34) 



Now, considering tlic functions Z, L and M as the reduced actions corresponding 
to the directions of the variables r, •& and (p, we can define a total reduced action 
of our system by summing the three reduced actions ^ 



5o(r, i9, p) = Z{r) + L(i?) + M(t?) . 
Taking up this equation in Eq. (34), we find 



(35) 



1 



2mo V 



So) - 



4too 



{5o,r} + l{5o,i?} + 



sin 1? 



{So,ip} 



+ V{r) 



8mo 8mo sin^ 



(36) 



Eq. (35) represents the QSHJE in 3-D for a central potential. At the classical 
limit {h^O), Eq. (35) goes to 



1 

2 m 



V,,^,^5o +V(r)=E. 



(37) 



which is the Classical Stationary Hamilton- Jacobi Equation (CSHJE). Remark 
that taking the classical limit in Eqs. (31), (32) and (33) make dL/d'd and 
dMjdtp) vanish, then, one cannot obtain the CSHJE (Eq.(37)). We deduce that 
the classical limit must be taken into Eq. (36). 

Note also that considering the 3D QSHJE Eq. (36), and separating variables 
with taking into accoimt Eq. (35), we get to Eqs. (31), (32) and (33) which 
lead to Eqs. (8), (9) and (10). After separating variables in Eq. (36), three 
constants of motion appear, they are the energy £, A and m;. 



The solutions of Eqs. (31), (32) and (33) are respectively 
Z{t) 



n arctan s ar — + br 
% 

h arctan •{ — -t- 6^ 



(38) 
(39) 



^ Which is also vahd for the classical case. 
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M {(fi) = harctani^a^Y > (40) 

Where a^, br, a^, b^}, a^p and a^p are real constants. Xi and X2 are two real 
independent solutions of Eq. (17). Ti and T2 are two real independent solutions 
of Eq. (18). Fi and F2 are two real independent solutions of Eq. (19). 

Now, knowing that: " the wronskian 

Ipl— ^^2-1— 

aq aq 

of two independent solutions of a second order differential equation with respect 
to g is a constant kg" , the expressions (38), (39) and (40) reduce to 

Z {r) — h a.Ycta.n I^Gr J + ' (^1) 
L(?9) = ?iarctan ja^; J ^ + b^^ , (42) 



M{(p) = ^arctan <{ / -^ + b^} , (43) 



dip 

where — krUr, — k^a^ and a^p ~ k^pUp are real constants. In Eqs. (41), 
(42) and (43), we need to determine only one of the two independent solutions 
of Eqs. (8), (9) and (10) respectively. 

Before entering into the elaboration of the quantum law of motion, we would 
like to exhibit the fact that for the central potential, only the radial function X 
depend on the form of the potential, otherwise the two angular functions T and 
F can be established without considering of the form of the central potential. 
Indeed, if we take a look on the table of solutions of the Schrodinger equation 
for such a potential we should write the angular function Ti as |10| 

^r' w - (1 - cos^(^))-""* [(1 - -omn , (44) 

and the angular function Fi as 

Fi((^) = sin(|m;|v5) if mi ^ 

= 1 if mi = . (45) 

It is usual to remark that we do not include the normalization constants in 
the two last expressions of the angular functions, since we consider the two 
mathematical solutions and not the physical's. 



3. The Quantum law of motion for a particle moving under the central 
potential 

Now let us investigate the quantum law of motion. In Ref. |7], we have 
presented a 3D formalism from which we described the quantum law of 3D 
motions for a general potential. We established a relation connecting between 
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the 3D velocity and the 3D quantum momentum (Eq. (5)). Using the spherical 
coordinates into Eq. (5) for the central potential case, we write 



.dZ 
dr 



d-d 



dM 

dip 



2[E-Vir)] 



(46) 



As in classic mechanics, we do not use directly Eq. (46) to plot the trajectories. 
We must first, separate the variables r, d and (p. This seems not possible from 
Eq. (46), and must be done in the 3D QSHJE (Eq. (36)). We noted in Sec. 
2 that the separation of variables in Eq. (36), after taking into account of Eq. 
(35), leads directly to Eqs. (31), (32) and (33). 

In Ref. 0, we have considered a deformation of the space like for general 
relativity. We treated the 3D quantum systems with a general potential by using 
the mechanics of a deformed space. This allowed us to introduce a quantum 
Lagrangian and derive the quantum law of motion (Eq. (5)). For the central 
potential, we proceed with the same manner. Thus, let us write Eq. (36) as 



/dzy 1 

\dr J arri^) 



dL 
dd 
1 



/dM 



r'^ si'n?{-d) \ dip J a^^{ip) 



1 



= 2mo[E - V{r)] , (47) 



where 



arr{r) 



1 - 



and 



1 



dZ 

dr 

dL\ 



/dMV 



{Z,r} 



(48) 



{M, ^} - 



are real functions depending on r, and ip respectively. They represents the 
non vanishing components of the metric tensor of the curved space created by 
the central potential. The Lagrangian is given as 



L{f, -d, ip, r, If) = — — arr{r) 



mo r 



mo sin^{'d) (p'^ 



(49) 



and leads to the dynamical equation 



mor^ mor^i?^ 

arr[rj H ai)i)[§) 



mo sin'^{'d) ip'^ 



2 w 2 ' 2 

After comparison between Eqs. (47) and (50), we get 

_ dZ/dr 
mor 



a^^{p) + V{r)=E . (50) 



(51) 
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dL/M 

= 7T > (52) 



mor^ sin^('!9)i^ 

Replacing these last equations into Eq. (47), we deduce directly Eq. (46). 
Now, replacing Eqs. (51), (52) and (53) into Eqs. (31), (32) and (33) respeo 
tively after taking into account of the relations (48), we find 



dZ_ . 

dr 



^ ' 2mor2 



(54) 



dL^^ l{l + l)n^ _ (mf - l/4)n^ ^ 
d'd mor^ mor^ sin^('i?) 



dM . (mj-l/m'^ 

^ = ^ — — 

dip mor^ sin (?9) 



The three above relations represent the decomposition of Eq. (46) with respect 
to the variables r, "& and <f. Indeed, if one takes the summation of Eqs. (54), 

(55) and (56) wc find Eq. (46). In fact, these equations describe the dynamical 
behaviour of a particle moving under the central potential, each one describes 
it in the corresponding coordinate (r, "& or (p). 

Note that only Eq. (54) depends on the nature of the central potential V{r) 
while Eqs. (55) and (56) do not depend on it. Then, Eqs. (55) and (56) can 
be studied whatever the nature of V(r). This is not the case for Eq. (54). So, 
as an application of the 3D quantum law of motion (Eqs. (46), (54), (55) and 

(56) ), we propose in the following sections to plot the trajectories of an electron 
moving under the action of the Hydrogen potential (Hydrogen atom). This will 
allow us to acquire a more comprehension of the 3D quantum law of motion. 

4. Classical trajectories of a particle moving under a central potential 

Before entering to plot the quantum trajectories, it is useful to present the 
classical Hamilton- Jacobi formulation of mechanics and also to establish the 
classical trajectories of a particle moving under the central potential. 

The Hamilton-Jacobi formulation is based on the CHJE given as 

i^^M\virl=E. (57) 

So represents the reduced action. Eq. (57) is written for the case of a central 
potential, after considering of the relation (35), as follows 



1 /dZy 1 1 fdlV 1 1 /dM^ ^ 



2mo \dr J ' 2mo \d^ J 2mo sin^ (^?) V df 
The separation of the variables into Eq. (58) gives 



dzV 



dr ) 



= 2mn 



E - V{r) - 



+ V{r) = E . (58) 



(59) 
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and 



• (61) 



a and /3 are real constants. As we know, the momenta dZ/dr, dL/d'd and 
dM/dip are connected to the speeds r, § and (f as follows 

— = mor , (62) 
^=morH, (63) 

dM 9 . 9 / „x . ,x 

— — = mor^ sm'^(i9)95 . (64) 
d(p 

Taking these equations into Eqs. (59), (60) and (61), we find 




r = Ml:^jE-Vir)-^, (65) 



1 / /32 

i cx — 

"mor^Y sin^(^) 



^ = ±—^Ja--r^ , (66) 



'^ = ± T^^' (67) 

TOor^ sm (77j 

The ± signs in Eqs. (65), (66) and (67) indicate that the motion of the particle 
may be in either direction on the r axis, and ip circles respectively. Now, 
we can integrate Eqs. (65), (66) and (67) to give the three time equation with 
respect to the coordinates r, ■& and (p, but before, we write these equations as 

/ 2 dr 
J—dt = ± , (68) 

dt . d'd 
2 = ± / = ' (69) 

dt dip 

- . (70) 



mor'^ sin (i9) (3 



We remark that to integrate Eq. (68) can be integrated directly after know- 
ing the form of the central potential (either analytically or numerically). How- 
ever, to integrate Eq. (69), we must first write r -already integrated from Eq. 
(68)- in function of the time t. Also, to integrate Eq. (70), we should set the 
variables r and d in function of t. In spite of the importance of the time equa- 
tions, we will focus on the spatial trajectory equations. These last's can be 
derived directly from Eqs. (68-70) after parameterizing the time t. Thus, we 
deduce 

sin^imip = ± , (71) 
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±- 



dr 



^a-/3Vsin2(,?) r^^/2moiE-Vir))-a/r 
dip = ± 



(72) 
(73) 



Integrating these last equations we find the trajectory equation !F{r, i?, ip) = 0. 

In order to compare the classical trajectories with the quantum trajectories, we 
plot in Fig. 1 the trajectory of an electron under the atomic potential (Hydrogen 
atom) 

nr-)=-— , (74) 

r 

where k = 9.10^5./ is the electric constant and e the charge of the electron. 
We give the energy E the value -(13.6/4) eV, which correspond in the quantum 
case to the energy of the first excited state of the Hydrogen atom (n = 2, (n is 
the principle integer)). We chose this value of energy in order to compare the 
classical trajectory (Fig. 1) with the quantum trajectories corresponding to this 
value of energy. For the same reason, we have chosen the constants a = 2 fi^ and 
(3 = (•\/3/2) h equal to the values of 1(1 + 1) and \/'mf — 1/4 h respectively, 
for Z = 1 and m; = 1. Then, the quantum trajectories corresponding to the 
classical trajectory are those plotted for n = 2, Z = 1 and mi = 1 (Figs. (8) and 
(10)). 

The classical trajectory showed in Fig. 1 is an inclined ellipse. This trajec- 
tory is obviously famous since it is the trajectory of a planet submitted to the 
gravitation force of the sun, and this force looks like the electric force of the 
atom. 



5. Quantum orbits of the Hydrogen's electron 

In this section we will plot the 3D quantum trajectories. First, let us write 
Eqs. (54), (55) and (56) as follows 



2 dt 



E - V{r) 



2mor^ 



dZ , 
— dr , 
dr 



dt ■■ 



{mf - 1/4) 



sin^(T9) 



dL 
dd 



dd 



mQr^s\v?{'d) dM , 

dt = — — — dip . 

{mf - l/4)h^ dip 



(75) 

(76) 
(77) 



Integrating Eqs. (75), (76) and (77), we find the time equations. In order to 
plot the quantum trajectories, we need to find the trajectories equations. In 
this order, we parameterize the time dt into Eqs. (75), (76) and (77) so as to 
obtain 



sin^(^) dM 
(mf - 1/4) ~d^ 



dip 



'^r\E-V{r))-l{l + l) 



dZ , 
— dr 
dr 



(78) 
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1(1 + 1) 



jmf - 1/4) 
sin^(i9) 



n -1 



2mo 



r2(£;-l^(r))-;(l + l) 



— dr 

dr 



dM 

dip 



dip 



1(1 + 1) 



K - 1/4) 



sin2(z9) - 1 



dL 

dd 



d§ 



(79) 



(80) 



In what follows Eqs. (78), (79) and (80) are using to plot the 3D quantum spatial 
trajectories (orbits). Expressions of the quantum momenta dZ/dr, dL/di) and 
dM /dip are deduced directly from the corresponding quantum reduced actions 
L and M (Eqs. (41), (42) and (43) respectively) 



dZ _^ 

dr (ttrXi J 



dL 

dd 

dM 

dip 



±- 



+ brXi)^ 



(a^i! f, + hai)^+Tl 



± 



ha<, 



(a^F, J^ + b^F,)^ . 



Ff 



(81) 



(82) 



(83) 



Knowing the expressions of the functions Xi , Ti and Fi , the quantum momenta 
will be known and the system of the three differential equations (Eqs. (78-80)) 
can be resolved numerically for the different values of the constants a^, for, 
h^, a^p and ht^. The ± sign in the expressions of the quantum momenta explicit 
the fact that the motion of the particle may be in either direction of the r axis, 
•d semicircle and ip circle. 

Note: In fact two of Eqs. (78-80) are sufficient to integrate and determine the 
quantum spatial trajectories since the third equation is automatically satisfied 
when the two others are. 

Now, let us consider the quantum motion of the Hydrogen's electron. We 
take up the expression of the atomic potential given in Eq. (74). In this case, 
the resolution of Eq. (14) is obvious and gives for the radial function X the 
following expression jlO| 



X['"'\r) 



exp 



r 
nar 



r 

an 



-I d 



7* n+l 



exp 



2r 
nan 



n being an integer number representing the principle quantum number, and 



(84) 



moke^ 



the Bohr's radius. The fimction X^ 



En 



(nl) 



corresponds to the energy 
2 fln 



(85) 



(86) 



Fixing the values of n, I and m into Eqs. (84), (44) and (45) and substituting 
these last equations into Eqs. (81), (82) and (83) respectively, then replacing 



12 



these lasts' and Eqs. (85) and (86) into Eqs. (78-80), we get a system of three 
differential equations which after integrating leads to the 3D quantum trajecto- 
ries of the electron for the corresponding bound state(n, /, mi). Naturally, these 
trajectories depend on the constants flr, 67-, 0^9, b^, a^p and htp. So, for several 
sets of those constants it correspond many trajectories (See Figs. (2-12)). In 
other words, for each bound state many trajectories exist exactly as we have 
found for the ID quantum motions. We will show this by plotting the trajec- 
tories of four bound states, the fundamental state (n = 1, Z = 0, m; = 0), and 
the excited bound states (n = 2, Z = 0, m; = 0), (n = 2, 1=1, mi = 0) and 
(n = 2, 1 = 1, mi = 1)^. With this aim, for each state we proceed as follows: 
First we explicit the fimctions T^^"^^ and Fi from Eqs. (84), (44) and 

(45). Secondly, we integrate numerically Eq. (75) and plot the time equation 
r{t) for different values of the constants 0^ and b^. After, by numerical methods 
we integrate two of Eqs. (78-79) for different values of the constants and plot 
the quantum spatial trajectories. 



a- the fundamental state 

For this state, we have n 
and Fi are deduced from Eqs. (84), (44) and (45) as 



For this state, we have n = 1, I = and mi = 0. The functions X^^^\ 7^*^°°^ 



Xr{r) = ^e.p(-^) , (87) 

Tl°°\d) = sini {§) , (88) 

F^"\^) = l. (89) 
Replacing Eqs. (87), (88) and (89) into Eqs. (81), (82) and (83), we find 

-2 ,,. 



^ = ± v_v (90) 

'^'^ (a,/^-'e^rfr + 6,)2 + l 

dL ^ _^ ha^ sin~^(^) 

rft? [atflog(|tan(|)|) + 6tf]2 + l ' ^ ' 



dip {a^ip + b^Y + 1 
Replacing Eq. (90) into Eq. (75), we find 



— dt = ± [-'-^ + —\ '!z_^!<l_:_ (93) 

^ """^ ^ ./(^) + +1 

Integrating numerically Eq. (93), we plot the radial time trajectories given in 
Fig. 2. In this figure, the electron moves between two radial extremities r = 
and r = 2 flo, and describes different trajectories for the same state. All these 

^Thc considered excited states in this article do not obeys to the selection law that we 
consider in spectroscopy of the Hydrogen atom. We just use these states in order to make a 
simple description of the quantum trajectories for different values of the quantum numbers n, 
I and m;. 
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trajectories pass throw some points constituting nodes with the same manner 
as for the ID motions [U 1^. The radial position of the nodes correspond 
to the zero of the radial function In particular, the nodes situated at 

r ~ 2ao correspond to the vanishing value of the radial velocity given by Eq. 
(93). Remark that all the trajectories oscillate about the trajectory for which 
the radial hidden variables take the values Ur = 0.36 and br = 0.52. 
Now, taking Eqs. (90), (91) and (92) into Eqs. (78), (79) and (80), we get 



4a<^ sin2(i?) , _ , "'■"°(~% + ^) 



—9 - 



^ ,dif:^± ^ 5 dr , (94) 



4a^ sin(i?) 



,log(|tan(f)|) 



■ dd = ±- ' ^ °^ — dr , (95) 



a. ^ a.sin-^W ^ 

(a^^ + 6^)2 + 1 ^ Klog(|tan(|)|)+6^]2 + l ^ ^ 

Integrating the last equations, we plot the 3D quantum spatial trajectories given 
in Fig. 3. This figure shows the 3D behaviour of the Hydrogen's electron. The 
electron describes different trajectories which pass from three points constituting 
spatial nodes. One of these nodes corresponds to r = and the two others to 
r = 2 Oq. 

It is useful to indicate that the motion of the electron in a state for which 
m; = is a purely quantum motion and cannot have a classical correspondent. 
The reason is because in Eqs. (60) and (61), /3'^ never takes a negative values 
while in Eq. (55) and (56), {mf — l/4)?i^ takes a negative values for to/ — 0, 
and this make impossible the construction of the classical corresponding of Eqs. 
(55) and (56) by given /J^ a negative values since /? is a real constant. 

b- excited state: n — 2, I — 0, mi — 

For this state, the radial function X^"^^^ is 

while the angular functions still given by Eqs. (88) and (89). Replacing Eq. 
(97) into Eq. (81), we find 



dZ 



-2 



6' 



d7=^7 ; , 72 ' (98) 

^^dr + bA +1 



^This position can be seen as the position of the nucleus in a first approximation and 
we may deduce that the electron fall on the nucleus. This is not the case because both 
electron and nucleus make a relative motion compared to their gravity center, so that the 

nucleus oscillate between the two radial positions r„uc = and r„uc = ^ ^j^j r^i , (for 

the hydrogen rnuc ~ 10"^'' X m). Then the electron has a real possibility to pass throw the 
position r,,; = without collision with the nucleus. 
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Taking Eq. (98) into Eq. (75), we find 

dt = ± -— + — ^ ^ ^ dr 



Too V a, 



(99) 

Integrating numerically Eq. (99), we plot the radial time equation (Fig. 4.) 
for this bound state. In this figure, the electron moves between two radial 

extremities r = and r = 8ao, and describes different trajectories for the same 
state. All these trajectories pass throw some points constituting nodes. The 
radial positions and 2ao of the nodes correspond to the zero of the radial 
function X^'^^K In particular, the nodes situated at r = Sao correspond to the 
vanishing value of the radial velocity given by Eq. (99). 
Taking Eqs. (98), (91) and (92) into Eqs. (78), (79) and (80), we get 



Aa^ sin^(z?) 
(a^ (fi + h^Y + 1 



4a^ sin(i?) 



dip = ± 2 I 

{^rl{^y\2-^y\^dr + br^ +1 

(100) 

= ± ; ° V <^o) aj — _ — 

^rl{^y\2-i;y'e^dr + br^ +1 

(101) 



,log(tan(f)) + M2 + l 



diD 1,9 sin ^ (■!?) 
7 r-^ d(p = ± — 5 — — d'& . (102) 

(a^<^ + 6^)2 + 1 ^ [a^log(tan(f)) + 6tf]2 + l ^ ^ 

Integrating Eqs. (100), (101) and (102) we plot the 3D quantum spatial trajec- 
tories given in Fig. 5 for this state. This figure shows clearly that for different 
values of the hidden variables the electron describes different trajectories which 
pass from three points constituting spatial nodes. One of these nodes corre- 
sponds to r = and the two others to r = Suo- Remark that the positions 
r = 2ao constituting nodes for the time radial trajectories do not constitute 
nodes for the spatial trajectories. We think that this is due to the angular be- 
haviour of the electron. In fact the angular trajectories in the {<p, i?) plan do 
not pass from any nodes for this state (see the Appendix). 

c- excited state: n = 2, Z = 1, mi = 

For this state , the radial function is 

Xl'^\^) = ^\^p(-^) , (103) 

do \ ^ O^o / 

and the angular function is deduced from Eq. (44) as 

r/^ °) (r) = sini (i?) cos(7?) , (104) 

without taking into account of the multiplicative constant (-2) obtained in 
the calculation of °^ since it has no influence on our equations. The angular 
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function still given by Eq. (89). Replacing Eqs. (103) and (104) into Eqs. 
(81) and (82), we find 



dZ 
dr 



dL_^ 



= ±- 



sin~^(t?) cos~^(t?) 



[a.(i±^+log(|tan(f)|))+5. 
Taking Eqs. (105) into Eqs. (75), we find 



(105) 



(106) 



^ ar 

— dt = ± 
mo 



+ 1 



dr 



e»o dr + br \ + 1 



(107) 



Integrating numerically Eq. (107), we plot the radial time equation (Fig. 6) 
for this bound state. The electron moves between two radial extremities r = 
(4 — 2 \/2) tto and r = {4 + 2 \/2) Aq, and describes different trajectories for the 
same state. All these trajectories pass throw some points constituting nodes. 
These nodes are situated at r = (4—2 \/2) ao and r = (4+2 \/2) ao corresponding 
to the vanishing values of the radial velocity given by Eq. (107). 
Taking Eqs. (105), (106) and (92) into Eqs. (78), (79) and (80), we get 



4a^ sin^(d) 
{a^ If + h^Y + 1 



d^p • 



± 



_2 



e »o rfr + 6r I + 1 



■ dr , 



(108) 



4atf (8 sin2(i?) + 1)-^ sin(i?) cos'^i^) 



+ cos(i?) 



+ log(|tan(t?/2)|) +h 



± 



■d-d = 



+ 1 



Co dr + br ] + 1 



dr , (109) 



(a^ip + b^)'^ + 1 



dip = ±' 



as (8 sin2(t^) + l)"i sin-i(i9) cos^^{^) 



(^7^+log(|tan(^/2)|))+6. 



■d'd . 



(110) 

Integrating Eqs. (108), (109) and (110) we plot the 3D quantum spatial trajecto- 
ries given in Fig. 7 for this state. This figure shows clearly that for different val- 
ues of the hidden variables the electron describes different trajectories which pass 
from two points constituting spatial nodes corresponding to r = (4 — 2 \/2) ao 
and r = (4 + 2 y/2) ao- 
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d- excited state: n = 2, Z = 1, mi = 1 

For this state , the radial function is given by Eq. (103), while the angular 
functions are 

T/")(t9) =sin5(i9) , (111) 

= sin(^) , (112) 

without taking into account of the multiplicative constant (-2) obtained in the 
calculation of t/"\ Replacing Eqs. (Ill) and (112) into Eqs. (82) and (83), 
we find 



dL 

dd 



= ±- 



ha^ sin (t?) 



dM 



(s^^+llog(|tan(f)|))+6. 



+ 1 



(113) 



(114) 



dip (a^ sm{ip) + cos{ip))'^ + cos((^)2 

The radial time equation is given by Fig. 6 since the radial function is the same 
for the present and the previous bound states. 

Taking Eqs. (105), (113) and (114) into Eqs. (78), (79) and (80), we get 



I aip siri?'{'d) 



{a^ sm{ip) + cos(<^))2 + cos2(i^) 



dip 



■ 



'dr + br] +1 



■dr, (115) 



(f sin^ii)) - 1) ^ sin"i(i9) 



■d-d : 



+ 1 



(s5^(fS^ + 5log(|tan(e/2)|))+&^' 



c^odr + br 



dr ,(116) 



■ dip = 



{a^ sm{ip) + 6<^cos(<^))2 + cos'^{ip) 

(I sin2(t?)-l)"^ sin-3(,?) 



±- 



(^5^ + ^log(|tan(f)|))+6.] +1 



■d-d 



(117) 



Integrating the last equations we plot the 3D quantum spatial trajectories. 
In Fig. 8, we compare a quantum trajectory integrated from Eqs. (115), (116) 
and (117) for the following values of the hidden variables o^. = 3.1, br = —0.4, 

= 0.6, 6tf = —0.2, = 1 and b^ = with the corresponding classical 
trajectory plotted from the classical equations (71), (72) and (73). We remark 
that the quantum trajectory seems to be, more or less, a deformation of the 
classical trajectory. To render this Note more obvious, we plot in Fig. 9 the 
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projection of the two trajectories into Y Z and XZ plans. The Fig. 9 shows 
that the quantum trajectory has a deformed eUiptic form. 

In Fig. 10, we plot different trajectories for different values of the hidden 
variables. These trajectories pass from two points constituting spatial nodes 
corresponding to r = (4 — 2 \/2) Oo and r = (4 + 2 \/2) Oo- 

6- Conclusion 

First, we would like to emphasize that this article is a continuation of the 
work already presented in Refs. 012 El B El El E], especially in Ref. 0. By this 
work, we want to bring a clarification and a checking of the results presented in 
Ref. [T]. Indeed, we have introduced in Ref. [7| a generalization to 3D of the 
deterministic approach of quantum mechanics presented in Refs. 0(2] ■ So, a 
3D quantum law of motion is established in Ref. [7] for a general form of the 
potential, and in the present article, we apply this law of motion for the central 
potential case and in particular to the Hydrogen's atom. We proceeded by three 
main steps, the first is the establishment of the 3D-QSIIJE and the quantum 
reduced action for the case of a central potential (Sec. 2). Then, in the second 
step, we established the 3D Quantum law of motion under the central potential 
(Sec. 3). Finally, in the third step, we plotted the 3D quantum trajectories of the 
Hydrogen's electron (Sec. 5). Analyzing our results, especially the trajectories 
obtained in Sec. 5, for different bound states of the electron, we find interesting 
outcomes. 

The first outcome is the fact that for the bound state case, the electron is 
trapped between two radial positions given by the two roots of the equation 

I mo 

and, if the electron reach an external position of the trapping zone delimited 
by the two roots, the electron will be ejected from the atom. This fact is well 
established in Fig. (11). Indeed, this figure shows that after passing the position 
r = 2 flo = 1.05834 x 10~^° m, which is one of the two roots of Eq. (118) for the 
bound state (n = 1, Z = 0, = 0) defining the trapping zone, the electron reach 
in a limited time a position of about lO^*' m. This result means that, after living 
the trapping zone the electron will be strongly ejected far from the atom. This 
is in disagreement with the standard quantum mechanics which suggest that for 
any position in the space the atomic electron has a non vanishing probability to 
be there. Our approach shows that the electron is trapped in a given trapping 
zone, otherwise it will be ejected far from the atom. Note that this result is 
even valid for all other excited bound states. A consequence of such a result is a 
partial localization of the electron in the space in the case of the bound states. 

The second outcome is the existence of nodes from which all the possible 
quantum trajectories plotted for the same state pass. In fact, we established 
the existence of nodes for two kind of trajectories, the radial time trajectories 
(Figs. (2), (4) and (6)) and the spatial trajectories (Figs. (3), (5), (7) and 
(10)). The radial time trajectories are similar to the ID time trajectories. So, 
in Ref. 0, we established that the nodes exist independently of the nature of 
the potential. This is also valid for the radial time trajectories and the nodes 
corresponds to the zeros of the functions A"*^"') and to the points where the 
radial velocity r vanishes. The spatial trajectories, in spite of their complexity. 
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pass throw points constituting nodes. The number and the positions of these 
nodes are estabhshed from Eq. (118). When the point r = is situated into 
the trapping zone defined by the two roots ri and r2 of Eq. (118), the number 
of nodes is three and their positions are r = 0, r — ri and r — r2- When the 
point r = is out of the trapping zone, the number of nodes is two and their 
positions are ri and r2. 

Another important outcome concerns the comparison between the quantum 
trajectories and the corresponding purely classical trajectory. We noted in Sec. 
5, that such comparison can be done only for the states for which m; 7^ 0, 
otherwise the classical conjugate momentum dM/d(p with respect to ip takes 
an imaginary value. This is the reason for what we take such a comparison 
only for the state for which n = 2, I — 1 and mi = 1. For this state, we find 
that either classical and quantum trajectories describe a closed curves in the 3D 
space which are contained between two radial positions r = (4 — 2 ^2) Oo and 
r = (4 + 2 \/2) Oq. The classical trajectory is an ellipse situated into an inclined 
plan (Figs. (1) and (8)), while the quantum trajectories have a complex forms 
and move in all the space. However, for some values of the hidden variables the 
quantum trajectories goes approximately to the classical trajectory with a an 
apparent residual deformation (Figs. (8), (9)). This is an obvious result since 
the quantum equations of motion reduce to the classical's. 

To conclude, we stress that our deterministic 3D approach of quantum me- 
chanics introduced in Ref. [3 and the establishment of the 3D quantum trajec- 
tories for a central potential (atomic potential) presented in the present article 
are a generalization of the deterministic ID approach of quantum mechanics 
introduced in Refs. CDIH- This deterministic approach, in the frame of its 
validity, will restore an objective reality, and permit a construction of a deter- 
ministic theory of quantum gravity. Such a construction will be an important 
step to start a unification of the fundamental interactions even the gravity. 

APPENDIX 

The angular behaviour of the Hydrogen's electron 

In this appendix, we plot many angular trajectories corresponding to several 
states of the Hydrogen's electron. (See Figs. (12)) 
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Figures 

Fig. 1: Classical 3D trajectory of the Hydrogen's electron. In order to consider 
the atomic scalars, we have chosen in Eqs (71-73), E = —13.6/4 eV, 
a — 2h^ and (3 — (\/3/2) . It is an ellipse, and we see clearly that the 
motion is done in a plan. 

Fig. 2: Quantum radial time trajectories for the Hydrogen's electron in the fun- 
damental state (E — —13.6 eV,7i = 1, ^ = 0, m; 0). For all the curves, 
we have chosen r{t = 0) = 2ao {do = 0.52917 . 10^^" m is the Bohr radius). 

Fig. 3: Quantum 3D spatial trajectories of the Hydrogen's electron in the funda- 
mental state {E = —13.6 eV, n = 1, / = 0, to/ = 0). 

Fig. 4: Quantum radial time trajectories for the Hydrogen's electron in the first 
excited state {E = —13.6/4 eV, n = 2,1 = 0, m; = 0). For all the curves, 
we have chosen r{t = 0) = 0. 

Fig. 5: Quantum 3D spatial trajectories of the Hydrogen's electron in the first 
excited state {E = -13.6/4 eY,n = 2,l = 0, mi = 0). 

Fig. 6: Quantum radial time trajectories for the Hydrogen's electron in the second 
excited state {E = —13.6/4 eV, n = 2,1 = l,mi = 0). For all the curves, 
we have chosen r{t = 0) = 0. 

Fig. 7: Quantum 3D spatial trajectories of the Hydrogen's electron in the second 
excited state {E = -13.6/4 eV, n = 2, I = 1, to/ = 0). 

Fig. 8: This figure shows a comparison between the purely 3D classical trajectory 
to the one of the possible 3D quantum trajectories of the Hydrogen's 
electron for the state {E = —13.6/4 eV, n = 2,1 = 1, mi = 1). For 
the corresponding values of the hidden variables (a^ = 3.1, fe,. = —0.4, 
= 0.6,5^ = —0.2, = 1, &^ = 0), we see easily from this picture, that 
the quantum trajectory can be considered as a deformation of the elliptic 
purely classical 3D trajectory.More comparison is done in Fig. 9, where 
projections of the two trajectories into YZ and XZ plans are done. 
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Fig. 9: Projection onto the YZ and XZ plans of the spatial trajectory for the state 

{E = -13.6/4 cV, n = 2, Z = 1, = \){ar = 3.1, K = -0.4,ffl,5 = 0.6, 
= —0.2, atp = l,bip = 0), and it corresponding purely classical trajec- 
tory, the quantum trajectory seems to be a deformation of the classical 
one. 

Fig. 10: Quantum 3D spatial trajectories of the Hydrogen's electron in the third 
excited state {E = -13.6/4 eV, n = 2, I = l,m; = 1). 

Fig. 11: The ejection of the electron from the atom when it live the trapping zone 
(0 < r < 2ao) in the case of the fundamental bound state {E = —13.6 
eV, n = 1, I = 0, mi = 0). We chose for the radial hidden variables the 
values Gr = 1.5 and br = —0.5. 

Fig. 12: The above figures show the relation between i? and cp in the electron tra- 
jectories for different / and mi. For many values of the hidden variables, 
we plot '!?(<p). (a): Z = 0, m; = 0. (b): Z = 1, m; = (c): Z = 1, m; = 1 . 
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